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ABSTRACT 

We report HETE-2 WXM/FREGATE observations of the X-ray flash, XRF 
020903. This event was extremely soft: the ratio log(S'x/5'^) = 0.7, where Sx and 

are the fluences in the 2-30 and 30-400 kcV energy bands, is the most extreme 
value observed so far by HETE-2. In addition, the spectrum has an observed 
peak energy E°^l^^ < 5.0 keV (99.7% probability upper limit) and no photons 
were detected above ~ 10 keV. The burst is shorter at higher energies, which is 

^Department of Physics, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8551, 
Japan 

^RJKEN (Institute of Physical and Chemical Research), 2-1 Hirosawa, Wako, Saitama 351-0198, Japan 
^Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, NM, 87545 
^Department of Astronomy and Astrophysics, University of Chicago, IL, 60637 

^Center for Space Research, Massachusetts Institute of Technology, 70 Vassar Street, Cambridge, MA, 
02139 

^Laboratoire d'Astrophysiquc, Observatoire Midi-Pyrenees, 14 Ave. E. Belin, 31400 Toulouse, Prance 

^Department of Physics, Aoyama Gakuin University, Chitosedai 6-16-1, Setagaya-ku, Tokyo 157-8572, 
Japan 

^National Astronomical Observatory, Osawa 2-21-1, Mitaka, Tokyo 181-8588, Japan 

^Tsukuba Space Center, National Space Development Agency of Japan, Tsukuba, Ibaraki, 305-8505, 
Japan 

^"Department of Astronomy, New Mexico State University, 1320 Prenger Mall, Las Cruces, NM, 88003- 
8001 

^^University of California at Berkeley, Space Sciences Laboratory, Berkeley, CA, 94720-7450 

^^Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, B.P.4346, 31028 Toulouse Cedex 4, Prance 

^^Department of Astronomy and Astrophysics, University of California at Santa Cruz, 477 Clark Kerr 
Hall, Santa Curz, CA 95064 

"NASA Goddard Space PUght Center, Greenbelt, MD, 20771 

^''Instituto Nacional de Pesquisas Espaciais, Avenida Dos Astronautas 1758, Sao Josedos Campos 12227- 
010, Brazil 

^^Department of Astronomy and Astrophysics, Tata Institute of Fundamental Research, Homi Bhabha 

Road, Mumbai, 400 005, India 

^'^Consigho Nazionale delle Ricerche (lASF), via Piero Gobetti, 101-40129 Bologna, Italy 
^^Faculty of engineering, Miyazaki University, Gakuen Kibanadai Nishi, Miyazaki 889-2192, Japan 



-3- 



similar to the behavior of long GRBs. We consider the possibihty that the burst 
hes at very high redshift and that the low value of -Epeak is due to the cosmological 
redshift, and show that this is very unlikely. We find that the properties of XRF 
020903 are consistent with the relation between the fluences S{7 — 30 keV) and 
S'(30 — 400 keV) found by Barraud et al. for GRBs and X-ray-rich GRBs, and are 
consistent with the extension by a decade of the hardness-intensity correlation 
(Mallozzi et al. 1995) found by the same authors. Assuming that XRF 020903 
lies at a redshift z — 0.25 as implied by the host galaxy of the candidate optical 
and radio afterglows of this burst, we find that the properties of XRF 020903 
are consistent with an extension by a factor ~ 300 of the relation between the 
isotropic-equivalent energy i^iso and the peak E'peak of the z/Fj, spectrum (in the 
source frame of the burst) found by Amati et al. for GRBs. The results presented 
in this paper therefore provide evidence that XRFs, X-ray-rich GRBs, and GRBs 
form a continuum and are a single phenomenon. The results also impose strong 
constraints on models of XRFs and X-ray-rich GRBs. 

Subject headings: Gamma rays: bursts (GRB 020903) 



1. Introduction 

Gamma-ray bursts (GRBs) that have a large fluence in the X-ray energy band (2-30 
keV) relative to the gamma-ray energy band (30-400 keV) are receiving increased attention. 
The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray 
Observatory detected 2704 GRBs (Paciesas et al. 1999). The spectra of 156 bright bursts 
exhibit a distribution of low-energy power-law indices a whose centroid is ~ —1, and a 
distribution of observed break energies E^^^^^ whose centroid is ~ 230 keV (Preece et al. 
2000), where E^^^^^ = (a - /5)(2 + a)~^E°ll^^. Here a, f3, and E°^l^ are the slope of the 
low-energy power-law index, the high-energy power-law index, and the energy of the peak of 
the pF,y spectrum of the Band function (Band et al. 1993), an expression that satisfactorily 
represents the spectra of almost all GRBs. In contrast, 36% of the bright bursts observed by 
GINGA have peak energies E^l^ in their photon number spectrum at a few keV and large 
X-ray to 7-ray fluence ratios (Strohmayer et al. 1998). 

The Beppo'&KX Wide Field Camera (WFC) detected events that are very similar to the 
soft GINGA GRBs; these events have been termed "X-ray Flashes" (XRFs) (Heise et al. 
2000).^^ The energy flux of these XRFs lies in the range 10~^ — 10~^ erg cm~^ s~^ and the 



Throughout this paper, we define "X-ray-rich" GRBs and XRFs as those events for which log[5x(2 — 
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low-energy photon index a of their spectra ranges from —3 to —1.2. The sky distribution 
of XRFs is consistent with isotropy, and there is no evidence that the sources are Galactic. 
The XRFs have tgo durations between 10 and 200 sec. The event rate of XRFs detected by 
the WFC is 3.4 events per year. Clarifying the connection between XRFs and GRBs could 
provide a breakthrough in our understanding of the prompt emission of GRBs. 

Kippen et al. (2002) made a detailed spectral comparison of GRBs and XRFs, using a 
sample of eighteen GRBs that were observed by BATSE and a sample of nine XRFs that were 
observed by both the WFC and BATSE. According to their joint analysis of WFC/BATSE 
spectral data, the low-energy and high-energy photon indices of XRFs arc —1 and ~ —2.5, 
respectively, which are no different from those of GRBs. On the other hand, XRFs have 
much lower values of E°^^ than do GRBs. Thus the only temporal or spectral difference 
between GRBs and XRFs appears to be that XRFs have lower E°^^ values. Kippen et al. 
therefore suggest that XRFs might represent an extension of the GRB population to events 
with low peak energies. Analyzing 35 HETE-2 GRBs seen by FREGATE, Barraud et al. 
(2003) demonstrate that the spectral properties of "X-ray rich" GRBs form a continuum 
with those of ordinary GRBs and suggest that XRFs may represent a further extension of 
this continuum. 

BATSE's low-energy threshold of ~ 20 kcV made it difficult for BATSE to detect XRFs. 
Ginga and BeppoSAX had the capability of detecting XRFs; however, Ginga could not 
determine the direction of the burst and the BeppoSAX GRBM had difficulty in triggering 
on XRFs. Consequently, these missions could not carry out in depth investigations of XRFs. 
In contrast, HETE-2 (Rickcr et al. 2003) has the ability to trigger on and localize XRFs, 
and to study their spectral properties, using the Wide-Field X-Ray Monitor [WXM; 2-25 
keV energy band; Kawai et al. (2003)] and the French Gamma Telescope [FREGATE; 6-400 
keV energy band; Atteia et al. (2003)], which have energy thresholds of a few keV. 

In this Letter, we report the detection and localization of XRF 020903 by HETE-2 
(Ricker et al. 2002) and present the results of a detailed study of its properties. Since this 
event was extremely soft and there was very little signal (a ~ 2a excess in the best selected 
energy range) in FREGATE, we focus our analysis on the WXM temporal and spectral data 
for the event. 



30 kev)/Sj{30 - 400 kev)] > -0.5 and 0.0, respectively. 
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2. Observations 
2.1. Localization 

XRF 020903 was detected with the HETE-2 WXM and the Soft X-ray Camera [SXC; 
0.5-10 keV energy band; Villasenor et al. (2003)] instruments at 10:05:37.96 UT on 2002 
September 3 (Ricker et al. 2002). The WXM flight locahzation was correct, but was not 
sent out because HETE-2 was pointing at the Galactic Bulge region at the time and WXM 
triggers were therefore not being sent to the GCN in order not to overwhelm the astronomical 
community with X-ray burst localizations. A GCN Notice reporting the localization of the 
burst, based on ground analysis (Graziani et al. 2003; Shirasaki et al. 2003) of the WXM 
data, was sent out 231 minutes after the burst. 

The WXM localization can be expressed as a 90% confidence circle that is 16.6' in radius 
and is centered at R.A. = 22M9"'25^ Dec. = -20°53'59" (J2000). A localization of the burst 
based on ground analysis (Monnelly et al. 2003) of the SXC data was distributed as a GCN 
Notice about 7 hours after the burst. Only a one-dimensional locahzation was possible using 
the SXC data, but this significantly reduced the area of the localization region for XRF 
020903. The improved localization produced by combining the SXC and WXM localizations 
can be described as a 90% confidence quadrilateral that is 4' in width and ~31' in length 
(see Figure 1). It is centered at R.A. = 22M9™01^ Dec. = -20°55'47" (J2000), and its four 
corners he at (R.A., Dec.) = (22M8'"48.00% -20°39'36.0"), (22M8'°33.60^ -20°42'36.0"), 
(22M9™10.80^ -21°10'12.0"), and (22M9"^30.00^ -21°10'48.0") (J2000). 

Detections of candidate optical and radio afterglows of XRF 020903, and the host galaxy 
of the candidate optical and radio afterglows, have been reported. Soderberg et al. (2002) 
discovered an optical transient within the HETE-2 SXC -|- WXM localization region at R.A. 
= 22*^48^42.34% Dec = -20°46'09.3", using the Palomar 200-inch telescope. These authors 
mention that the optical transient brightened by ~ 0.3 — 0.4 magnitudes between about 7 and 
24 days after the XRF, and suggest tliat the re-brightening might be due to an associated 
supernova. However, the optical transient apparently faded by over a magnitude only three 
days later (Covino et al. 2002). Spectroscopic observations of the optical transient, using the 
Magellan 6.5m Baade and Clay telescopes, detected narrow emission lines from an underlying 
galaxy at a redshift z = 0.25 ± 0.01, suggesting that the host galaxy of the optical transient 
is a star-forming galaxy [Soderberg et al. (2002); see also Chornock & Filippenko (2002)]. A 
fading bright radio source at the position of optical transient was detected using the Very 
Large Array (Berger et al. 2002). Hubble Space telescope observations of the XRF 020903 
field reveal the optical transient and show that its host galaxy is an irregular galaxy, possibly 
with four interacting components (Levan et al. 2002). These detections likely represent the 
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first discoveries of the optical and radio afterglows, together with the host galaxy, of an XRF. 

In our analysis of the prompt emission of XRF 020903, we apply a "cut" to the WXM 
photon time- and energy-tagged data (TAG data), using only the photons from the pixels 
on the three wires in the X-detectors (XAO, XAl, and XA2) and the four wires in the Y- 
detectors (YAl, YA2; YBO, YBl) that were illumnated by the burst and that maximize the 
S/N of the burst light curve, in the same manner as we did for GRB 020531 (Lamb et al. 
2003). We use this optimized TAG data when performing our temporal and spectral analyses 
of this event. 



2.2. Temporal Properties 

Figure 2 shows the light curve of XRF 020903 in four WXM energy bands. The time 
history of the burst in the 2-5 and 5-10 keV energy bands has two peaks. Clearly, there is 
no significant flux above 10 keV. Table 1 gives the tso and tgo durations of the burst in the 
2-5 keV, 5-10 keV, and 2-10 keV energy bands. The duration of the burst is longer in the 
lower energy band; this trend is similar to that seen in long bright GRBs (Fenimore et al. 
1995). 

2.3. Spectrum 

As we have seen, the light curve of XRF 020903 shows two peaks: the first occurring 
in the time interval 0-8 s, and the second occurring in the time interval 8-13 s. The S/N of 
the first peak is much higher than that of the second. In addition, inspection of the burst 
light curve in the 2-5 and 5-10 keV energy bands suggests that the second peak is much 
softer than the first. For these reasons, we analyze the spectrum of the burst in three time 
intervals: 0-8 s, 8-13 s, and the total duration of the burst, 0-13 s. The background region 
we use is 40 seconds in duration and starts 45 seconds before the burst. 

The WXM detector response matrix has been well-cahbrated using observations of the 
Crab nebula (Shirasaki et al. 2002). In the spectral fits, we include only the photons that 
registered on the four wires in the X-detectors and the five wires in the Y-detectors that 
were illuminated by the burst, as mentioned above. Since the variation in the gain is not 
uniform at the ends of the wires in the WXM detectors (Shirasaki et al. 2000), we use only 
the photon counts that registered in the central ±50 mm region of the wires to construct 
the spectra of the burst. We include all of the photons that register in the central regions of 
these wires (i.e, we use the full 2-25 keV energy range of the WXM). The relation between 
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pulse height and energy in the WXM is non-hnear and is different for each wire. In order 

to extract the strongest possible constraints on the parameters of the spectral models we 
consider, we treat each individual WXM wire separately but take the normalizations on all 
wires to be the same. For the same reason, we do not re-bin any of the pulse height channels 
in the WXM and in the FREGATE, and we carry out a set of fits for the total duration of 
the burst (0-13 s) that include the spectral data from both the WXM and the FREGATE. 
We use the XSPEC vll.2.0 software package to do the spectral fits. 

Table 2 presents the results of our time-resolved and time-integrated spectral analysis of 
the burst. In this analysis, we consider the following models: (1) blackbody, (2) power-law 
model, (3) power-law times exponential (PLE) [the COMP model in Preece et al. (2000)], 
and (4) Band function (Band et al. 1993). The Galactic value of Nh in the direction of the 
burst is 2.3 x 10^° cm~^ (Dickey & Lockman 1990), which is negligible (i.e., it is undetectable 
in the WXM energy range). Furthermore, the WXM and FREGATE data do not request 
Nh as a free parameter (i.e., introducing Nh as an additional free parameter produces only a 
small change in x^)- There is therefore no need to include Nh as a parameter in the spectral 
fits, and we do not do so. In the Band model fits, we have fixed a = —1, which is a typical 
value for GRBs, in order to better constrain the remaining parameters. All of the models 
provide acceptable fits to the data; i.e., the data do not request models more complicated 
than a blackbody or a power-law. 

However, essentially all GRB spectra arc well- described by the Band function (Band et 
al. 1993; Preece et al. 2000), and the analysis of Kippcn et al. (2002) shows that at least 
some XRF spectra are also well-described by the Band function. Fits to the WXM data for 
all three time intervals using the power- law model give spectral slopes a < —2 with high 
significance. For example, comparing the minimum value xLin — ^3.4 corresponding to the 
best-fit value of the spectral slope a — —2.8 and the value — 60.6 at a = — 2 for the 
power-law fit to the average spectrum of the burst (i.e., the time interval 0-13 s), we find that 
Ax^ = 7.2 for one additional parameter. Thus, using the Maximum Likelihood Ratio Test, 
a < —2 at the 99.3% confidence level. From this evidence, we conclude that the peak i?pcak 
in vFi, lies near or below 2 keV, the lower limit of the energy range of the WXM detectors. 

There is evidence of spectral softening between the first and second time intervals. In 
particular, a power-law fit to the first time interval gives ai = — 2.4to;6 ^^^^ Xmin,i = 75.1, 
while a power-law fit to the second time interval gives a2 = —4. 2^37 and Xmin2 — 71.1. In 
contrast, a power-law fit to the first and second time intervals, but with a = ai = a2 gives 
a = — 2.86]!;q;82 and xLin — 152.2. The first (more complicated) model includes the second 
model as a special case (i.e., the models are nested). Comparing xLin the two models, 
we find that Ax^ = 6.0 for one additional parameter. Thus, using the Maximum Likelihood 
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Ratio Test , there is evidence of spectral softening at the 98.6% confidence level. 

Figure 3 shows a comparison of the observed count spectrum and the count spectrum 
predicted by the best-fit power-law model, for the time intervals 0-8 s and 8- 13s. Figure 4 
shows the same comparison, except for the total duration of the burst (0-13 s). Table 3 gives 
the peak photon number and energy fluxes (in 1 s) and the fluence of XRF 020903, assuming 
the power-law model. 

Using the power-law model that best fits the burst-averaged WXM plus FREGATE 
spectral data, we find fiuences Sx{2 - 30 keV) = 8.2^2 j x 10^^ erg cm-^ and 5^(30 - 
400 keV) = l-6l!li'3 ^ 10~^ erg cm~^ where the quoted errors give the 90% confidence intervals. 
Thus the ratio of fluenccs \og[Sx{^ - 30 keV)/^^(30 - 400 kcV)] = 0.7, with a 90% lower 
limit of 0.3, making this burst not only an XRF but the most extreme example of an XRF 
observed so far by HETE-2. 

A comparison of the power-law and Band function fits to the first peak, which has a 
much higher S/N than the second peak, provides modest evidence for an ^'pglj^ near 2 keV, the 
lower limit of the energy range of the WXM detectors. Specifically, we find thatAx^ = 4.34 
for one additional parameter, which means that the data requests the (more complicated) 
Band function model at the 89% confidence level. However, the evidence is clearly not of 
high statistical significance, and in this fit we fixed a at —1, its typical value for GRBs. 

We therefore choose to place an upper limit on E^l^. The appropriate model to use is 
the Band function, since (as we have already mentioned) the spectra of almost all GRBs and 
at least some XRFs are well-described by this function. However, this presents a problem: 
the Band function has two distinct ways of representing a power-law spectrum in the detector 
energy range. First, it can do so by having E'break — > 0, so that only the high-energy, pure 
power-law part of the Band function is visible in the energy range of the detector. Second, 
it can do so by having E'break oo and £'o ^ oo, where Eq is the "cutoff energy" of the 
cutoff power-law that constitutes the low-energy part of the Band function. In this limit, 
the limiting power-law is actually the cutoff power-law, but the cutoff energy is so large that 
the curvature of the model is imperceptible in the detector energy range. 



"We compute the peak photon number flux in the WXM 2-5, 5-10, and 2-10 keV energy bands, using 
the best-fit power-law model parameters for the average photon energy spectrum of the burst, and the ratio 
2.731 of the photon flux in the 1 s time interval containing the largest number of photons and the average 
photon flux in the - 13 s time interval. Wc compute the peak photon energy flux in the WXM 2-5, 5-10, 
and 2-10 keV energy bands in exactly the same way, except that we use the ratio 3.247 of the total photon 
energy flux (found by weighting each photon with its energy and summing the energies) found in the 1 s 
time interval containing the largest total photon energy and the average photon energy flux in the - 13 s 
time interval. 
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We solve this problem by developing a new statistical method. This method uses a 
constrained Band function which is parameterized by two quantities, and p. The 

constrained Band function is perfectly able to make both pure power-law spectra and power- 
law times exponential spectra of the required curvature in the detector energy range, but only 
the high-energy part of the Band function is allowed to produce a pure power-law spectrum. 
We describe this new method in detail in Appendix A. This method has general applicability 
to all instruments when the spectra of the bursts considered have E°^ly. near or below the 
low-energy threshold of the detector. It is necessary to demonstrate that the photon index 
P < —2 before applying this method. 

In applying the constrained Band function method to XRF 020903, we jointly fit the 
WXM and the FREGATE data. Figure 5 shows the posterior probability density distri- 
bution for ^'peak ^li^t we find using this approach. From this posterior probability density 
distribution, we find a best-fit value = 2.7 keV, that 1.1 keV < E°^^ < 3.6 keV 

with 68% probabilty, and that < 4-1 and 5.0 keV with 95% and 99.7% probabihties, 

respectively. 

We conclude that the properties of XRF 020903 are very similar to those of long GRBs, 
with the exception that the observed peak energy -Speak ~ 3 keV is ~ 100 times smaller. 
The extremely low value of E°^^ seen in XRF 020903 is similar to the smallest value found 
among the 9 XRFs whose spectra were determined by jointly fitting BeppoSAX WFC and 
BATSE data (Kippen et al. 2002). 



3. Discussion 

3.1. Source Properties 

We exclude the possibility that XRF 020903 is a Type 1 X-ray burst (XRB) on the 
following grounds. First, its galactic latitude is 6 = —61.5° (using the center of the combined 
WXM plus SXC error box) , and there is no known persistent X-ray source or globular cluster 
in this error box. Since Type I X-ray burst sources lie in the Galactic plane or in globular 
clusters, and have persistent X-ray emission, XRF 020903 is unhkely to be an X-ray burst 
on locational grounds alone. Second, the time history of XRF 020903 is not FRED-hke (i.e., 
it does not exhibit a fast rise and an exponential decay), while those of XRBs typically are. 
Third, although the blackbody model gives an acceptable fit to the spectra of the first and 
second peaks in the time history of XRF 020903, the derived blackbody temperatures are ~ 
1.0 keV. These temperatures are lower than those of almost all Type 1 X-ray bursts [which 
typically have temperatures T f» 2 keV; see, e.g. Lewin, van Paradijs & Taam (1993)]. For 
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these reasons, we conclude that XRF 020903 is an XRF and not an XRB. 

The extremely low value of E°^l^ observed for XRF 020903 is remarkable. If the observed 
spectrum of XRF 020903 were the redshifted spectrum of a typical GRB, the implied redshift 
would be z ~ 100, using the best-fit value of E°^l^ = 2.7 keV observed for XRF 020903 and 
the mean value of ^'brlak ^^^r the sample of 5500 spectra formed from the brightest 156 BATSE 
GRBs (Preece et al. 2000). A redshift of this magnitude would be hard to understand, and is 
certainly not expected if long GRBs are associated with the collapse of massive stars (Lamb 
& Reichart 2000). It is also wildly inconsistent with the measured redshift z = 0.25 of the 
host galaxy (Soderberg et al. 2002) of the candidate optical (Sodcrbcrg ct al. 2002) and radio 
(Berger et al. 2002) afterglows of XRF 020903. It is therefore difficult to attribute the low 
observed value of -Bpelk XRF 020903 to cosmological redshift. 



3.2. Fluence and Pecik; Energy Correlations 

In Figure 6, we plot XRF 020903 in the (S'3o-4oo55'7-3o)-plane, where Sjso and 5*30-400 
arc the energy fiuences of the bursts in the 7-30 and 30-400 keV energy bands. For the 
value of 5*7-30 and 5*30-400, we use the best-fit power- law model for the average spectrum of 
WXM and FREGATE. Also plotted in this figure are the 35 GRBs whose spectra have been 
determined using HETE-2 FREGATE data (Barraud et al. 2003). Figure 6 shows that the 
properties of XRF 020903 are consistent with the relation between 5*7-30 and 5*30-400 found 
by Barraud et al. (2003). 

In Figure 7, we plot XRF 020903 in the (52-4oo,^peak)-plane, where E°^l^ is the peak of 
the observed uF,, spectrum. For -Bpeak, we plot the 99.7% upper limit (5.0 keV). The proper- 
ties of XRF 020903 are consistent with an extension by two decades of the hardness-intensity 
correlation (Mallozzi et al. 1995; Lloyd-Ronning, Petrosian & Mallozzi 2000) between 5*30-400 
and found by Barraud et al. (2003). 

Amati et al. (2002) demonstrated that there is a relation between Eiso and the burst- 
averaged value of £^peak (i-e., £^peak for the time-averaged spectrum of the burst). Assuming 
that the candidate optical and radio afterglows of XRF 020903 are indeed the afterglows of 
XRF 020903, and therefore that the redshift of the underlying host galaxy is the redshift of 
the XRF, we can calculate the isotropic-equivalent radiated energy E[so and the upper limit 
on the burst-averaged peak energy -Epeak of the uF^, spectrum in the source frame, in the 
same way as did Amati et al. (2002). Figure 8 shows that the properties of XRF 020903 are 
consistent with an extension by a factor of ~ 300 in Eiso of the relation found by Amati et 
al. (2002). 
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Figures 6-8 provide evidence that XRFs, "X-ray-rich GRBs," and GRBs form a con- 
tinuum, and are therefore the same phenomenon. 



3.3. Constraints on Theoretical Models of XRFs 

A variety of theoretical models have been proposed to explain XRFs [see, e.g., Zhang 
& Meszaros (2003) for a comparative discussion of several of these models] . In the off-axis 
GRB jet model (Yamazaki, loka, & Nakamura 2002, 2003), XRFs are the result of viewing 
the jet of an ordinary GRB off-axis, so that relativistic beaming shifts the 7-rays into the 
X-ray range. In the clean fireball model, XRFs are due to the relativistic pair plasma in 
the GRB jet becoming optically thin much later than usual, at which time the relativistic 
bulk Lorentz factor T has already decreased to a relatively low value (Mochkovitch et al. 
2003). In the dirty fireball model, XRFs occur when there is significant baryon loading of 
the GRB jet, so that F never reaches large values (Dermer, Chiang, & Bottcher 1999; Huang, 
Dai, & Lu 2002). In the universal jet model, XRFs are the result of viewing the GRB jet 
off-axis, where F is lower because of the structure of the jet (Rossi, Lazzati, & Rees 2002; 
Woosley, Zhang, & Heger 2003; Zhang & Meszaros 2002; Meszaros, Ramirez-Ruiz, Rees, 
& Zhang 2002). In the uniform jet model, the different properties of XRFs, "X-ray-rich" 
GRBs, and GRBs are due primarily to different jet opening angles, with larger jet opening 
angles associated with lower values of F (Lamb, Donaghy & Graziani 2003). 

Any such model of XRFs must reproduce the correlation found by Barraud et al. (2003) 
between Sj^^o and 5*30-400) and the evidence we report in this paper for correlations between 
5'2-4oo and E^l^, and especially, £^iso and £^peak ~ the latter spanning nearly five decades in 



4. Conclusions 

In this paper, we have reported HETE-2 WXM/FREGATE observations of the X-ray 
fiash, XRF 020903. This event was extremely soft: the spectrum had a best-fit peak energy 
-^pcak ~ 2.7 keV and -Ep^^k < ^.0 keV (99.7% probability upper limit) and no photons were 
detected above ~ 10 keV. The burst is shorter at higher energies, which is typical of long 
GRBs. We considered the possibility that the burst lies at very high redshift and that the 
low value of E°ll^ is therefore due to the cosmological redshift, and showed that this is 
very unlikely. We find that the properties of XRF 020903 are consistent with the relation 
between 5'7_3o and 5'3o_4oo found by Barraud et al. (2003) for GRBs and X-ray-rich GRBs, 
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and are consistent with an extension by two decades of the hardness-intensity correlation 
(Mallozzi et al. 1995; Lloyd-Ronning, Petrosian & Mallozzi 2000) between 5*30-400 and E'^^l-^ 
demonstrated by the same authors. Assuming that XRF 020903 hes at a redshift z = 0.25 
as imphed by the host galaxy of the candidate optical afterglow of this burst, wc find that 
the the properties of XRF 020903 are consistent with an extension by a factor ~ 300 of the 
relation between Eiso and i^peak in the source frame of the burst found by Amati et al. (2002) 
for GRBs. When combined with earlier results, the results reported in this paper provide 
strong evidence that XRFs, X-ray-rich GRBs, and GRBs form a continuum and are a single 
phenomenon. The correlation found by Barraud et al. (2003) between S^sq and <5'3o_4oo, and 
the evidence we find in this paper for correlations between 5'3o_4oo and E^^^, and especially, 
Eiso and -Bpeak, provide strong constraints on any model of XRFs and X-ray-rich GRBs. 
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A. The "Constrained" Band Function For Soft GRBs 

In the spectral analysis of GRBs, one occasionally encounters events (such as XRF 
020903) that are so soft that they present themselves as pure power-laws with power-law 
index /3 < — 2 in the energy range of the detector. The natural interpretation of such spectra 
is that the break energy S'break separating the two functional parts of the Band function is 
near or below the lower boundary of the detector energy range. 

This situation creates a problem for fits of the Band function, in that the Band function 
has two distinct ways of conforming to a power-law in the detector energy range: 

1- ^break 0, SO that Only the high-energy, pure power-law part of the Band function is 
visible in the energy range of the detector. 

2. -Ebrcak ~^ oo, Eq oo, whcrc Eq is the "cutoff energy" of the cutoff power-law that 
constitutes the low-energy part of the Band function. In this limit, the limiting power- 
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law is actually the cutoff power-law, but the cutoff energy is so large that the curvature 
of the model is imperceptible in the detector energy range. 

Therefore, despite the fact that the numerical value of the power-law index is such that 
we are certain that we should be dealing with the high-energy part of the Band function (i.e., 
the index is < —2), the low-energy part of the function can "horn in" on the fit, altering the 
physical inferences drawn from the spectrum. 

This situation is particularly a problem for the estimation of -Epeak- Since we know that 
we are in case 1, we also know that we ought to have at least a firm upper limit on E^cak, 
since -Epeak is always necessarily less than inbreak; which is at the low end of the detector 
energy range. On the other hand, the case 2 limit implies -Epcak oo. Unfortunately, the 
data don't care which side of the Band function makes the power-law, so no discrimination is 
possible between the two cases. Consequently, we can't constrain ii^peak at all using a normal 
Band function fit. 

The approach we have chosen to deal with this situation in the case of XRF 020903 
is to fit a constrained Band function to the data. That is, we consider a three-dimensional 
subspace of the full four-dimensional Band function parameter space, choosing the subspace 
with a view to satisfying the following criteria: 

1. It is perfectly possible to make both pure power-laws and cutoff power- laws of the any 
desired curvature in the detector energy range. 

2. Only the high-energy part of the Band function is allowed to produce a pure power-law. 

We define the three-dimensional subspace in the following way: consider a Band function 
parametrized by low- and high-energy indices a and f3, and by a cutoff energy Eq. The well- 
known relation between Eq and -Ebreak is -E'break = (a — P)Eo. We impose the constraint 
condition on our family of fitting functions 

-E^break = -E^pivot X (-EoZ-Epivot) ^, (Al) 

where £^pivot is some suitably chosen energy, in the general neighborhood where the GRB has 
appreciable emission. E'break and Eq are then inversely related, and are equal to each other 
when both are equal to £^pivot- 

When Eq < i^pivot, then ii^break > -Epivot, and the function is essentially a cutoff power- law 
in the energy range of interest. 
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On the other hand, when Eq > E'pivot, then E'break < -Epivot, and as £'0 — > oo, £'break 0. 
In other words, when the low-energy part of the Band function is trying to imitate a power 
law, the break energy becomes small enough to force the low-energy part of the function 
below the energy range of interest, where it cannot be seen and therefore can do no harm. 
Any pure power-law work must thus be done by the high-energy part of the Band function. 

The resulting spectral function has three parameters (including the scale), rather than 
four. The two input shape parameters can be chosen arbitrarily from the set {a, f3, Eq, inbreak, 
-E'peak}- The remaining parameters may then be determined by algebraic relationships. 

We have found it most convenient to adopt Epeak and /3 as our parameters. The choice 
of -E'peak is dictated by the necessity of estimating its value, or at least an upper bound 
on its value. The choice of f3 is convenient because one may then impose the parameter 
bound f3 < —2, which guarantees that the formal expression for -Epcak niay be meaningfully 
interpreted as the energy of the peak of the i^Ei, distribution. This bound on f3 is an important 
part of the specification of the fitting family of models. Were it not imposed, it would be 
possible for the formal expression for Ep^ak to exceed E'break, so that at large values of £^peak 
the fit could always produce a /3 > — 2 power-law in the detector energy range. The result 
would be an extended tail of constant for arbitrarily large values of E'peak- 

Figure 9 shows the constrained Band function, with /3 = —2.5 and -Epivot = 4 keV, 
for different values of E'peak- This figure shows that £'peak increases, S'break also necessarily 
increases, so that Eq is forced to smaller and smaller values by the constraint, which increases 
the curvature (and the value of a) . 

Figure 10 shows the constrained Band function, with -Epeak = 4 keV and -Bpivot — 4 keV, 
for different values of /3. The progression from some curvature at low energy (/3 = —2.0) 
to almost none (/3 = —4.0) is evident, as is the fact that as the curvature disappears, the 
resulting power-law is produced by the high-energy part of the Band function. 

Figures 9 and 10 show that the constrained Band function is perfectly able to make both 
pure power-laws, and cutoff power-laws with any desired curvature in the detector energy 
range. Figure 10 demonstrates that, in the constrained Band function, a power-law spectrum 
is always produced by the high-energy part of the Band function. 

The choice of -Bpivot is dictated by the following considerations: 

1. Epiyot must be low enough to prevent the low-energy part of the Band function from 
making a power- law in the energy range of interest. If £'pivot were 1 GeV (say), then 
the Band function would have no difficulty making Eq large and a < — 2, which is 
what we are trying to prevent by introducing the constraint. So £'pivot should be "as 
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low as possible." 

2. £^pivot must not be so low that we cannot adequately fit any curvature that may exist 
in the spectrum. If E'pivot were 1 eV (say), then whenever E'break was in or above the 
energy range where the spectrum is appreciable, Eq would be so tiny that the curvature 
of the model would be huge, much too large to fit the data well. 

One way of choosing £^pivot is to calculate its value using the best-fit parameters from 
a fit of a free Band function, using -Epivot = (.E'o-E'break)''^'^^- This choice, which effectively 
chooses the unique constrained subspace of the full parameter space that contains the best- 
fit free Band function, allows the constrained family of functions to optimally fit whatever 
curvature the data may seem to hint is required. 

We must require that the inferences that we draw from the spectral fit should be robust, 
in the sense that they should not depend strongly on the specific choice of fi'pivot- So the 
proper use of this constrained Band model involves not only choosing a representative value 
of -Epivotj but also varying Epi^ot in some reasonable range, to make sure that the conclusions 
about parameter estimates and bounds are unaffected by the choice of -Epivot- 

Figure 11 shows the constrained Band function, with Epcak = 4 keV and (3 = 2.5, for 
different values of -Epivot- Once again, as the low-energy curvature disappears, the resulting 
power-law is produced by the high-energy part of the Band function. Figure 11 also shows 
that the shape of the spectrum in the detector energy range is insensitive to the specific 
choice of £'pivot, within a reasonable range. Thus the conclusions about parameter estimates 
and bounds are unaffected by the choice of fi'pivot- 

Figure 12 shows the constrained Band functions with parameters that best fit the 13 s 
spectrum of XRF 020903, for different fixed values of i^pivot- This figure illustrates the fact 
that the shape of the best-fit model is essentially unchanged in the energy range of the WXM 
for choices of E'pivot within a reasonable range.' 

Finally, we give the algebraic relationships necessary to recover the remaining Band 
function parameters assuming that -Epeak and P are given. Let x = £'peak/-E'pivot- Then, 



In Equation (A2), we have resolved the ambiguity in the choice of root of a quadratic equation 
by requiring that when P + 2 < 0, then a -|- 2 > 0, so that -Bpeak is in fact the peak energy 
of the uFi, distribution. 



a 




(A2) 
(A3) 



Eo = (2 + a)^peaic. 
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Table 1: Temporal properties of XRF 020903. 



Energy Band 






(keV) 


(s) 


(s) 


2 - 5 


5.8±0.9 


10.6±0.2 


5 - 10 


2.4±0.2 


4.3±2.2 


2 - 10 


4.9±0.6 


9.8±0.6 



Note. — The quoted errors correspond to ±la. 



Table 2. Results of fits to the spectrum of XRF 020903. 



Time region 


Model 


kT 
(keV) 


a 






rpohs 
peak 

(keV) 


Xi (DOF) 


0.0-8.0 


blackbody 
power-law 
cutoff power-law 


-'^■"-'^-0.20 


_2 4+0-5 
-1.0 (fixed) 






q 1+1.9 
"^•-■--1.1 


1.08 (62) 
1.21 (62) 
1.14 (62) 




Band 




-1.0 (fixed) 


< 


-2.4 


o 4+1.7 


1.16 (61) 


8.0-13.0 


blackbody 
power-law 
cutoff power-law 




-4 2+1-^ 
-1.0 (fixed) 






< 2.0 


1.13 (62) 
1.15 (62) 

1.14 (62) 




Band 




-1.0 (fixed) 


< 


-3.3 


< 2.0 


1.15 (61) 


0.0-13.0 


blackbody 
power-law 
cutoff power-law 


0.871°:?° 


-1.0 (fixed) 






2 4+^-2 


0.79 (62) 
0.86 (62) 
0.81 (62) 




Band 




-1.0 (fixed) 


< 


-2.7 


2 4+^-2 


0.82 (61) 


0.0-13.0 


blackbody 


0.901°:?^ 










0.85 (177) 


(with FREGATE) 


power-law 
cutoff power-law 




-2 6+°-'' 

^•"-0.5 

-1.0 (fixed) 








0.86 (177) 
0.85 (177) 




Band 




-1.0 (fixed) 


< 


-2.3 


< 4.1 


0.86 (176) 



Note. — The quoted errors correspond to the 90% confidence region 
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Table 3: Peak photon number and energy fluxes (in 1 s) and fluences in various energy bands 
for XRF 020903. 







2-5 keV 


5-10 keV 


2-10 keV 


Peak flux (ph cm"^ s""^) 
Peak flux (10~^ ergs cm~^ s' 
Total fluence (10~^ ergs cm" 




1.9 ±0.7 
10.411? 
4.2 ±0.9 


S3+°-^^ 

1 7+0.8 
-•-■'-0.7 


2.2 ±0.8 
14.7 ±5.3 
5.9 ± 1.4 



Note. — All of the quantities in this table are derived assuming a power-law model 
for the spectrum. The quoted errors correspond to the 90% confidence region. 
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Fig. 1.— The HETE-2 WXM/SXC localization for XRF 020903. The circle is the 90% 
confidence region for the WXM localization and the belt-like region is the portion of the 
90% confidence region for the one-dimensional SXC localization that lies within the WXM 
90% confidence circle. The final localization is the intersection of the WXM and SXC 
localizations (Ricker et al. 2002a). The point labeled "OT" is the location of the candidate 
optical (Soderberg et al. 2002) and radio (Berger et al. 2002) afterglows of XRF 020903. 
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Fig. 2.— The light curve of XRF 020903 in four WXM energy bands: 2-5 keV, 5-10 keV, 
10-25 keV, and 2-25 keV (top to bottom). The hght curve is binned in one second bins. The 
vertical dotted lines show the 0-8 and 8-13 second time intervals bracketing the first and 
second peaks of the burst light curve. We have performed model fits to the spectra of the 
burst during these two time intervals, and to the entire duration of the burst (0 -13 seconds). 
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Fig. 3. — Comparison of the WXM spectra for the time intervals 0-8 and 8-13 s. The 
observed (crosses) and predicted (histogram) count rates are shown in a different color for 
each of the nine WXM wires that we have included in the fits. The spectral model is a power 
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Fig. 4. — The WXM and FREGATE spectra for the entire time interval 0-13 s. The observed 
(crosses) and predicted (histogram) count rates are shown in a different color for each of the 
nine WXM wires that we have included in the fits. The spectral model is a power law with 
fixed photoelectric absorption (see Table 2). 
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Fig. 5. — The posterior probability density distribution for -Epeak- The vertical solid lines 
define the 68% probabihty interval for -Epeak^ while the dashed and dotted lines show the 
95% and 99.7% probability upper limits on E°ll^. We find a best-fit value E^]^l^ = 2.7 keV, 
that 1.1 keV < E^^l^ < 3.6 keV with 68% probability, and that E^^l^ < 4.1 and 5.0 keV 
with 95% and 99.7% probability. 
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Fig. 6. — The (5'7_3o,5'3o_4oo)-plane, showing the location of XRF 020903, using the total 
fluence in the 7-30 keV energy band and in the 30-400 keV energy band (filled circle). The 
crosses are the locations of the 35 HETE/FREGATE GRBs studied by Barraud et al. (2003). 
The sohd fine is the relation, 87^30 = 3 x lO'^ S|]o%o> found by Barraud et al. (2003). The 
properties of XRF 020903 are consistent with this relation. 
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Fig. 7.— Tlie (52-4oo,^peak)-plane, showing the location of XRF 020903. For E°^l^, we plot 
the 99.7% probability upper limit (5.0 keV). We calculate the 2-400 keV fluence using the 
best power-law model fit jointly to the WXM and the FREGATE data. The crosses show 
the locations of 12 of the HETE-2 GRBs studied by Barraud et al. (2003) for which E°^l^ is 
relatively well determined. The properties of XRF 020903 are consistent with an extension 
by two decades of the hardness- intensity correlation found by Barraud et al. (2003). 
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Fig. 8. — The (£'iso,-E'peak)-plane, where ii^iso is the isotropic-equivalent radiated energy be- 
tween 1-10^ keV and -Epeak is the peak of the uF^, spectrum, both measured in the rest frame 
of the burst. The filled circle in the lower left-hand corner is the location of XRF 020903. 
The ten open circles are the BeppoSAX GRBs reported by Amati et al. (2002). The solid 
line is given by the equation, Epeak = 89{Ei^o/l0^^erg)^-^ keV. The properties of XRF 020903 
are consistent with an extension of this relation by a factor ~ 300. 
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Fig. 9. — Constrained Band functions, with (3 = —2.5 and ii^pivot = 4 keV, for different 
values of -Epeak- All functions have been normalized to 1 keV~^ at 10 keV. The two vertical 
hues at 2 keV and at 25 keV show the WXM bandpass. The spectra shown are (decreasing 
monotonically from the top at low energy), -Bpeak =1 keV, 2 keV, 4 keV, 6 keV, and 8 keV, 
respectively. As E^cak increases, £^break also necessarily increases, so that Eq is forced to 
smaller and smaller values by the constraint, increasing the curvature and the value of a. 



-30- 




E (keV) 

Fig. 10. — Constrained Band functions, with Epeak — 4 keV and E'pivot = 4 keV, for different 

values of f3. All functions have been normalized to 1 keV~^ at 10 keV. The two vertical lines 
at 2 kcV and at 25 keV show the WXM bandpass. The spectra shown are for j3 = -2.0, -2.5, 
-3.0, -3.5, and -4.0, which can be distinguished by the increasing steepness of their slopes 
at high energy. The progression from some curvature at low energy {f3 = —4.0) to almost 
none {(3 = —2.0) is evident, as is the fact that as the curvature disappears, the resulting 
power-law is produced by the high-energy part of the Band function. 
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Fig. 11. — Constrained Band functions, with -Epcak = 4 kcV and f3 = —2.5, for different 
values of -Epivot- All functions have been normalized to 1 keV^^ at 10 keV. The two vertical 
lines at 2 kcV and at 25 keV show the WXM bandpass. The spectra shown are (increasing 
monotonically at low energy) for -Epivot = 2 keV, 4 keV, 6 keV, and 8 keV, respectively. Once 
again, as the low-energy curvature disappears, the resulting power-law is produced by the 
high-energy part of the Band function. Note also that the shape of the constrained Band 
function is insensitive to the specific choice of E'pivot within a reasonable range. 
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Fig. 12. — Constrained Band functions with parameters that best fit the 13 s spectrum 
of XRF020903, for different fixed values of -Epivot- The two vertical lines at 2 keV and at 
25 keV show the WXM bandpass. All functions have been normahzed so that the integral 
from 2 keV to 25 keV is one photon. The five spectra shown in the plot corresponding to 
-Spivot = 4 keV, 5 keV, 6 keV, 7 keV, and 8 keV (the 7 keV and 8 keV largely overlap each 
other). This figure illustrates a robust aspect of the constraint procedure: the best-fit model 
is essentially unchanged in the WXM spectral band despite a factor-of-two change in the 
value of -Bpivot- 



